Second-order perturbations of a zero-pressure cosmological medium: 
comoving vs. synchronous gauge 
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Except for the presence of gravitational wave source term, the relativistic perturbation equations of a 
zero-pressure irrotational fluid in a flat Friedmann world model coincide exactly with the Newtonian 
ones to the second order in perturbations. Such a relativistic-Newtonian correspondence is available 
in a special gauge condition (the comoving gauge) in which all the variables are equivalently gauge 
invariant. In this work we compare our results with the ones in the synchronous gauge which has 
been used often in the literature. Although the final equations look simpler in the synchronous 
gauge, the variables have remnant gauge modes. Except for the presence of the gauge mode for the 
perturbed order variables, however, the equations in the synchronous gauge are gauge invariant and 
can be exactly identified as the Newtonian hydrodynamic equations in the Lagrangian frame. In 
this regard, the relativistic equations to the second order in the comoving gauge are the same as the 
Newtonian hydrodynamic equations in the Eulerian frame. We resolve several issues related to the 
two gauge conditions often to fully nonlinear orders in perturbations. 



I. INTRODUCTION 

The general relativistic cosmological linear perturba- 
tion theory was first developed by Lifshitz in 1946 [1]. 
Lifshitz took the synchronous gauge condition in which 
the perturbations of the time-time part and the space- 
time part of the metric tensor are equal to zeros; this 
gauge condition can be taken to fully nonlinear order 
without losing any physical degree of freedom [2]. The 
synchronous gauge condition has been popular in the cos- 
mological perturbation literature despite the complicat- 
ing fact that, except for the zero-pressure case, there are 
remnant gauge modes for both the spatial and temporal 
gauge conditions. There exist other spatial and tempo- 
ral gauge conditions which fix the gauge transformation 
property completely in general situation, thus without 
any remaining gauge mode [3-5] . This point was clarified 
by Bardeen [6,7]. In a zero-pressure medium the density 
perturbation equation in the synchronous gauge coincides 
with the one in the comoving gauge [1,5]. The density 
perturbation equation in the comoving gauge condition is 
known to resemble the Newtonian equation most closely 
[5,6], and the equations coincide in the zero-pressure case 
[1,8]. Thus, in the zero-pressure case the density pertur- 
bation equation in the synchronous gauge coincides with 
the Newtonian one to the linear order [1,8]. 

The synchronous gauge was also used in the nonlin- 
ear perturbation studies [9], and Kasai [10] has derived 
second-order differential equations for density perturba- 
tion which is valid to fully nonlinear order. Although, 
such an equation in the synchronous gauge naturally has 
proper linear limit which corresponds to the Newtonian 
equation, it has been unclear whether such a correspon- 
dence continues to the nonlinear situation. Recently, we 
have successfully shown an exact relativistic-Newtonian 



correspondence of scalar-type perturbations to the sec- 
ond order based on the comoving gauge [11-13]. In the 
zero-pressure case our comoving gauge condition differs 
from the conventional synchronous gauge in the spatial 
gauge condition. In this work we will investigate the case 
in the original synchronous gauge. We will show that al- 
though the equations in the synchronous gauge look sim- 
pler than the ones in the comoving gauge, the variables 
still have remaining (spurious) spatial gauge mode to the 
second order. The equations in the synchronous gauge, 
however, are gauge invariant and can be identified as the 
Newtonian hydrodynamic equations in the Lagrangian 
frame. Whereas, the equations in the comoving gauge 
can be identified as the Newtonian hydrodynamic equa- 
tions in the Eulerian frame. 

Results in Sees. II and the Appendices are valid to fully 
nonlinear order in perturbations, and unless mentioned 
otherwise results in the remaining sections are valid to 
the second order in perturbations. We closely follow no- 
tations used in [11-13]. We set c = 1. 



II. FULLY NONLINEAR PERTURBATIONS 

The energy conservation equation and the Raychaud- 
hury equation give [14,12,13] 



fj, + fi0 = 0, 

9 + \e 2 + & ab & ab 



UJ ah UJab + ^Gfl 



A = 0, 



(1) 

(2) 



where 6 = u a . a is an expansion scalar based on a fluid 
four- vector u a ; a a b and uj a b are the shear and the rota- 
tion tensors based on u a , respectively; tildes indicate the 
covariant quantities and the Latin indices indicate space- 
time components. An overdot with tilde is a covariant 
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derivative along the u a vector, e.g., A = A,aU a - By com- 
bining these equations we have 

/ r \ 7 / r \ 2 

I A | 1 / M | ~afc~ , ,~.ab,~, 



A / 3 \ // 



a a6 a a6 + u ab u ab - 4ttGA + A = 0. 

(3) 



Equations (l)-(3) are fully nonlinear and covariant, thus 
valid to all orders in perturbations. 



A. Temporal Comoving Gauge 

In this work we will assume an irrotational fluid, thus 
u) a b = 0. We will consider two different gauge conditions. 
In both gauge conditions we will have 



u a = 0, 



(4) 



due to a common temporal gauge condition together 
with the irrotational condition; the Greek indices indi- 
cate space components. If we introduce the spatial part 
of the four-vector as 



(5) 



where i>a^ is a vector-type perturbation (thus trans- 
verse), the irrotational condition sets Va^ = and our 
temporal comoving gauge sets v = 0. Since u a = the 
fluid four-vector in this gauge coincides with the normal 
frame four- vector h a with h a = 0. Notice that our tem- 
poral comoving gauge condition v = (together with the 
irrotational condition) implies u a = 0. This differs from 
the ordinarily known comoving frame condition which 
sets u a = [15]. In our case the normalized (u a u a = —1) 
fluid four-vector u a becomes 
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M = 



r,00 ■ 



u a = 0; 



u° 



VH™, u a = - 



~ g0a 



(6) 



In the zero-pressure case the momentum conservation 
equation implies g ao = —1/a 2 where a is the cosmic scale 
factor of the Friedmann background world model. In the 
ADM approach [16], our temporal comoving gauge it = 
together with the irrotational condition implies vanish- 
ing momentum vector J a = —h^T^ = 0. The ADM mo- 
mentum conservation equation in Eq. (13) of [11] gives 
ff„ = where g 00 = -l/N 2 , thus N = N(t). In another 
way, since the acceleration vector a a = u a -bu b = (In N)^ a 
vanishes (i.e., geodesic flow) for the zero-pressure irrota- 
tional flow, we have N = N(t); see Eqs. (27) and (42) 
of [11]. Without losing generality we can set N = a(t). 
Thus we have 



1 

7^- 



(7) 



Thus, Eq. (6) becomes 

u = —a, ii a = 0; u° = -, u a = —ag 0a , (8) 

which is valid to fully nonlinear order. We can show 
that to all orders in perturbations the fluid quantities are 
independent of the spatial gauge condition which could 
affect g 0a ; see the Appendix A. 



B. Nonlinear perturbed equations 

We introduce perturbations 

jl = li + 5fjt, 9 = 3H-k, (9) 

where H = a/a and S = Sfi/n; an overdot denotes a time 
derivative based on background proper-time t. The 9 is 
an expansion scalar of the fluid four-vector which is the 
same as the normal four-vector because u a = in our 
case. Using Eq. (8) we have 



A = A (1 + 5) + fi ( S - -N a S, 
a 



6 = 3H 



—N a n , 



(10) 



where N a is the shift vector in the ADM notation with 
JV a = a 2 g 0a ; the spatial indices of the ADM variables 
are based on h a p = g a p. 
Equations (1), (2) give 



(A + SHfi) (1 + S) 



5--6, a N a -(1 + S)k 



= 0, 



(11) 



3 [H + # 2 ) + 4ttG/i - A 



k - — K a N a + 2Hk - iTiGfiS - \n 2 - d ab d ab 



= 0. 

The background parts give 

A + 3i/^ = o, z(h + h 2 

The perturbed parts give 
1=(1 + S)k, 



k + 2Hk = -k 2 + a ah a ah + AirC^S, 



(12) 



AirGfi - A = 0. (13) 



(14) 
(15) 



where 5 = 5 — a 1 (5. ct A^ a . By combining these equations 
we have 



S + 2HS - 47rG/i<5 = 4ttG/z<5 



2 , 4 (S) 2 



31 + 5 



(l + 5)~a ab ~a ab . 

(16) 
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These equations are valid to the fully nonlinear orders 
in perturbations, subject only to the temporal comov- 
ing gauge condition, the zero-pressure condition, and the 
irrotational condition. 



C. The synchronous gauge 

Under the synchronous gauge we set goa = (thus 
N a = 0) using the spatial gauge condition (together with 
the irrotational condition), thus 



li = fi = (i. 
Thus, Eqs. (14)-(16) simply give 
5={1 + 6)k, 



(17) 

(18) 
(19) 



k + 2Hk =\n 2 + o ah b ah + AnGfiS, 
S + 2HS- 4itGh5 = 4^G/i<5 2 + + (1 + 6) ~o ah ~o ah , 

6 1 + 

(20) 

which are valid to the fully nonlinear order. Using A = 
5/(1 + S) Kasai [10] has derived 

2 A 2 

A + 2HA - AnGfiA — 1 " ^~ ab ~ 



3 1 - A 



(1 - A) b a0 ~o ah . 



(21) 



To nonlinear order in perturbations the above equations 
are incomplete yet because of a ab a a b term. Later we will 
show that these equations in the synchronous gauge dif- 
fers from the equations in the comoving gauge to the sec- 
ond order. Furthermore, although these equations look 
simple, we will show that 5 (thus A as well) and k still 
have remnant gauge modes to the second order. In Sec. 
Ill B we will show that to the second order the equations 
are gauge invariant and can be identified with the New- 
tonian hydrodynamic equations in the Lagrangian frame. 
In this regard, the equations in the comoving gauge corre- 
spond to the Newtonian hydrodynamic equations in the 
Eulerian frame. 



III. SECOND-ORDER PERTURBATIONS 



(3) 

order variables are based on g a p, and a vertical bar in- 
dicates the covariant derivative based on g^ ; g^ could 
become 5 a p in a flat Friedmann background. We ignored 
the transverse vector- type perturbation variables. We in- 
troduce x = a(/3 + a-y). The perturbed variables can be 
regarded as nonlinearly perturbed ones to any order in 
perturbations. 

To the second order, from Eqs. (55), (57) of [11] we 
have 



N a = -(3' a , 
~o ah ~o ah = KpK@ = 



(23) 



where N a is evaluated to the linear order; is a trace- 
free part of the extrinsic curvature. We note that a ah a a b 
is spatially gauge invariant to the second order, see Sec. 
IV. 

Before comparing equations in the two different spatial 
gauges, we compare our v in Eq. (5) with the notation 
used in [11] to the second order in perturbations. In [11] 
we introduced the fluid quantities based on the normal- 
frame vector h a and provided the relation of fluid quanti- 
ties between the energy-frame (E) and the normal-frame 
(TV). Our fluid four- vector u a is based on the energy 
frame which sets q a = 0. The energy-frame fluid four- 
vector is introduced in Eq. (53) of [11], and using the 
relations given in Eqs. (87), (88) of [11] we have 



B a + AB a 



* = a (U Q £ 

= J «L L_ 

\fi+p [n+p) 2 



(<J/i + 5p) QZ + Q p N n af} 



(24) 



Using the decomposition of the normal-frame flux vector 
Q% = {ji+p){-v. a +v ( a ] ) in Eq. (175) of [11] and setting 
= we have 

v,a = v, a ^— [(Sfi + 5p) v. a + u ,/3 n a/3 ] . (25) 

H + p 

Thus, the temporal comoving gauge v = in [11] implies 
v = and vice versa. 



As the metric we take 



ds 2 = -a 2 (1 + 2a) drj 2 - 2a 2 /3^dr]dx a 
+a 2 \gV> (1 + 2<p) + 2 1MP + 2C% 



dx a dx p , (22) 



where a, (3, 7 and ip are spacetime dependent perturbed- 
order variables, and C^p is a transverse and tracefree 
perturbed-order variable. Spatial indices of perturbed 



A. The comoving gauge 

In [11-13] we took the temporal comoving gauge and 
the spatial 7 = gauge 



v = 0, 7 = 0. 



(26) 



In this work, we call this the comoving gauge. Thus, we 
have /3 = xl a - 



3 



The momentum conservation equation in Eq. (105) of 
[11] gives 



a 



2a 2 



(27) 



Thus, apparently, a does not vanish to the second order. 
Later we will show that if we take (3 = as the spatial 
gauge condition instead of 7 = 0, we have vanishing a. 
However, we prefer 7 = as the spatial gauge condi- 
tion because it fixes the spatial gauge degree of freedom 
completely (as long as we simultaneously take the tem- 
poral gauge which removes the temporal gauge degree of 
freedom completely, like our v — 0), see Sec. VI of [11]. 
Whereas, (3 = fails to fix the spatial gauge degree of 
freedom completely, thus having remaining gauge degree 
of freedom even after imposing the gauge condition, see 
Sec. IV B. 

In our gauge the fluid four-vector in Eq. (8) becomes 



u = -a, 

_n 1 



u a = 0; 

r ^ = \r 



(1 - lip) 5% - 2C 



(t)a 



(28) 



Thus, although we prefer to call this the temporal comov- 
ing gauge (see [6,7]), because u a = and u a 7^ 0, our 
fluid four-vector corresponds to the normal four-vector 
rather than the comoving one. 

Using Eq. (23), Eqs. (14), (15) give 



K = Sk, 



a/3 



(29) 



(30) 



These also follow from the energy-conservation equation 
and the trace part of ADM propagation equation in Eqs. 
(104), (102) of [11]. 

In [12,13] we identified to the second order 



Sfi = Sg, k = — — V • u, 



(31) 



where 5g and u are Newtonian density and velocity per- 
turbations, respectively. As we ignore the rotational 
mode, the velocity is of potential type with u = Vit. 
Apparently, we need \ to the linear order only, and to 
that order we have [12,13] 



(32) 



where we assume a flat Friedmann background world 
model. With these identifications of the relativistic met- 
ric and energy-momentum perturbation variables (these 
are equivalcntly gauge-invariant combinations, see Sec. 
IV B) with the Newtonian hydrodynamic variables, Eqs. 
(29), (30) give 



5 + -V • u = --V • (Su) , 

a a 

-V • (u + Hu) + AitGliS = -^V • (u • Vu) 



(33) 



(34) 



By combining these we have 



5 + 2-5- AirGfiS = [oV • {5u)\ + -Ly • (u • Vu) 



a 2 dt 



(35) 



which also follows from Eq. (16). Except for the presence 
of the gravitational waves as source terms Eqs. (33)- (35) 
are valid exactly in the Newtonian system [17]. 

Although our relativistic equations are valid to the sec- 
ond order, the Newtonian equations are valid to fully 
nonlinear order. Thus, all nonvanishing higher-order per- 
turbation terms in the relativistic case are pure general 
relativistic corrections. Recently, we have presented such 
pure general relativistic correction terms appearing in the 
third order perturbations in [18]. 



B. The synchronous gauge 

The synchronous and comoving gauge conditions cor- 
respond to taking [2] 



0, (3 = 0; a = 0. 



(36) 



In this work, we call this simply the synchronous gauge. 
Thus, we have 7 = x/a 2 . If we take v = and (3 = 
as the temporal and the spatial gauge conditions, respec- 
tively, the momentum conservation equation gives a = 
to all orders in perturbations; although this is well known 
in [2], we give proofs in the Appendix B. Thus, in the 
zero-pressure medium without rotation we can simultane- 
ously impose the comoving (v = 0) and the synchronous 
(a = 0) temporal gauge conditions as long as we also 
take (3 = as the spatial gauge condition [2] ; Kasai took 
these conditions in his work in [10]. 

The original synchronous gauge used by Lifshitz [1] 
took a = and (3 = as the temporal and the spatial 
gauge conditions, respectively. These gauge conditions 
are known to be incomplete in fixing both the tempo- 
ral and the spatial gauge modes even to the linear or- 
der. Thus, even after imposing these gauge conditions 
we have remaining gauge modes present in the solutions, 
in general. Meanwhile, v = and 7 = fix the tem- 
poral and spatial gauge degree of freedoms completely, 
thus no gauge mode is present in the solution, see Sec. 
IV. Since the original synchronous gauge implies v = 
(the nonvanishing solution of v is the remnant tempo- 
ral gauge mode) in the zero-pressure case, we only have 
to pay attention to the possible presence of the spatial 
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gauge mode. In this gauge we have <?oo — —a 2 = l/<? 00 
and g a = = g 0a . Thus, the fluid four- vector in Eq. (8) 
becomes 



u 



-a, u c 



0; u° = -, u a = 0, 
a 



(37) 



which can be compared with Eq. (28) in the comoving 
gauge. Thus, since u a = 0, our fluid four-vector corre- 
sponds to the conventionally known comoving four- vector 
[15], and simultaneously normal because u a — as well. 
All the statements in the above two paragraphs are valid 
for all perturbational orders. 

Using Eq. (23), Eqs. (18), (19) give 



5 — k = 5k, 

k + 2Hk - 4-rrGfiS 



72 X, a p + C, 



a/3 



(38) 



(39) 



These also follow from the energy-conservation equation 
and the trace part of ADM propagation equation in Eqs. 
(104), (102) of [11]. By combining these equations we 
have 



6 + 2H5 - AnGfiS = 5 2 + AttG^S 2 



a/3 



(40) 



Apparently, these equations in the synchronous gauge 
look simpler than Eqs. (29), (30) in the comoving gauge. 
Compared with Eqs. (29), (30) in the comoving gauge, 
in Eqs. (38), (39) we lack the convective-derivative-like 
terms in the left-hand-sides. By changing the time 
derivatives as 



9 d 1 ,„ , „ d 1 „ 



dt dt 



(41) 



we can show that Eqs. (38)- (40) become the same ones 
in the comoving gauge in Eqs. (29), (30), and (35); in the 
last step of Eq. (41) we used Eq. (32) which is valid for 
a flat background. 

If we make the same identification of the density and 
velocity perturbations as in Eqs. (31), (32), thus assum- 
ing a flat background, Eqs. (38), (39) become: 



1 1 
5 + -V • u = — SV ■ u, 

a a 

-V • (u + Hu) + AnG^S = ~ (V'V) (V/ju a ) 



(42) 



(43) 



Ignoring the gravitational waves, these equations can be 
identified as the Newtonian hydrodynamic equations in 
the Lagrangian frame. 

Although equations in the synchronous gauge look sim- 
pler than the ones in the comoving gauge, the presence 



of additional convective-like terms in the comoving gauge 
allows us to make exact (except for the gravitational 
waves) correspondence with the Newtonian hydrody- 
namic equations in the Eulerian frame [12,13]. Whereas, 
the equations in the synchronous gauge can be identi- 
fied as the Newtonian equations in the Lagrangian frame. 
However, the variables in the synchronous gauge still 
have the remnant spatial gauge mode due to incomplete 
fixing nature of the spatial gauge condition [3 = in that 
gauge. That is, to the second order, 5 and k in the syn- 
chronous gauge have the remaining gauge modes, see Sec. 
IV B. 

Now, we can relate the variables in the synchronous 
(S) gauge to the ones in the comoving (C) gauge. From 
Eqs. (59), (56), and (32) we have 

Ss = S C +(^J ^Xdt + V75,Gaugc^) • VJc, 

K S =K C + (^J ^Xdt + V75,Gaugo^ • Vk c , (44) 

where 7s,Gau g e is the gauge mode present to the linear 
order in 7; see the next section. In a flat background, 
from Eq. (32) we have V% = au. Notice that, even after 
ignoring the gauge modes 5s and ks naturally differ from 
5c and kc, respectively, because the final equations are 
different. Using Eq. (44), Eqs. (38)-(40) give Eqs. (29), 
(30), and (35). 

Although the variables in the synchronous gauge have 
remnant spatial gauge mode, somehow the equations 
in the synchronous gauge coincide with the Newtonian 
ones in the Lagrangian frame. Meanwhile, the Newto- 
nian hydrodynamic equations have nothing to do with 
the gauge mode which appears only in the relativistic 
treatment. We can show that the situation is consis- 
tent in the synchronous gauge. From Eqs. (54), (55) the 
gauge mode of 5s. Gauge = £ a V Q <5c is proportional to the 
linear-order solution of 5c; similarly, the gauge mode of 
K 5, Gauge = £ a V a Kc is proportional to the linear-order 
solution of Kc- Thus, the behaviours of the gauge mode 
cannot be distinguished from the solutions to the lin- 
ear order, and can be absorbed to the linear order solu- 
tions. We can also check that the gauge modes in Eq. 
(44) cancel out in Eqs. (38) and (39). In this sense, Eqs. 
(38) and (39), thus Eqs. (40), (42) and (43) as well, are 
gauge- invariant. Therefore, to the second order in the 
synchronous gauge, although the variables have remnant 
gauge mode, the equations are gauge invariant; this hap- 
pens because the gauge mode temporally behaves exactly 
like one of the physical solutions. 

A similar situation occurs to the linear order in the 
original synchronous gauge which took only a — = (3 
[1]. Under these gauge conditions Lifshitz derived 



5 + 2H5 - AttG^S = 0, 



(45) 



which is the LHS of Eq. (40) and concides with the later 
derived Newtonian equation [8]. However, under these 
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gauge conditions (i.e., without taking v = 0), 5 still have 
the remnant gauge mode due to the incomplete fixing na- 
ture of the temporal gauge condition a = 0. It happens 
that the temporal behaviour of gauge mode of 6 is propor- 
tional to H which coincides with one of the two physical 
solutions, see [19]. Thus, although 5 has the gauge mode 
Eq. (45) is gauge invariant. In our synchronous gauge 
which also takes v = the temporal gauge condition is 
fixed completely, but a similar situation repeats due to 
an incomplete fixing nature of the spatial synchronous 
gauge condition (/3 = 0) now to the second order in per- 
turbation. 



IV. GAUGE ISSUE 
A. Gauge transformation 

Under a transformation between two coordinates x a — 
x a + £ a , the gauge transformation properties of all metric 
and energy-momentum variables to the second order are 
presented in Sec. VI of [11]. 

Since both the synchronous gauge and the comoving 
gauge take w = 0we have 

i° = 0. (46) 

This follows from Eqs. (234) or (238) of [11]: in the nor- 
mal frame by setting Q a = (i.e., v = 0) in both gauges 
we have £° a = 0; or in the energy frame, by setting 
V a - B a +AB a + 2V l3 C afj = (i.e., v = 0) in both 
gauges, we again have £° a = 0. Thus, without losing 

generality we can take £° = 0. 

To the second order, with |° = 0, from Eqs. (229), 
(232) of [11] we have 

a = a - a a C - P, a ^' - \e"&, 

S = 6-S^ a , K = K-K, a £ a , (47) 

where £ Q = £ a with the index of £ a based on ; a prime 
indicates the time derivative based on the conformal time 
n with dr/ = dx° = dt/a. The gauge transformation prop- 
erty of k follows from the scalar nature of the expansion 
scalar 8 with 9 = 3H — k where 9 is based on the normal 
frame; for the gauge transformation of a scalar quantity, 
sec Eq. (239) of [11]. To the linear order, from Eq. (252) 
of [11] we have 

0, a = 0, a +&, 7,a = 7,a - &« • (48) 

Thus, x = a (/?+7') is gauge invariant to the linear order, 
and 

a - a a r a + l(3, a 0' a , 5-5, a y a , K-K, a r a , (49) 
a-a aT a - (p+l^'J l[ a , (50) 
are gauge invariant to the second order. 



B. Two gauges 

In the comoving gauge, by imposing 7 = in all coor- 
dinates (i.e., 7 = = 7), from Eq. (48) we have 

= 0. (51) 

Thus, the spatial gauge transformation property is fixed 
completely. From Eq. (47) we have 

a = a, 5 = 6, k = k, (52) 

and each variable in this gauge has unique gauge- 
invariant counterpart as S and k in Eq. (49) and a in 
Eq. (50). Thus, we can equivalently regard all variables 
in this gauge as (spatially and temporally) gauge invari- 
ant ones. For example, <5„ j7 = S — 8. a j' a is a unique 
gauge invariant combination which is the same as 6 in 
the v = = 7 gauge conditions; for an explicit form 
of £„ >7 including v, see Eq. (282) in [11]. We note that 
these results (i.e., values remain the same in the comov- 
ing gauge conditions, complete fixing of the gauge de- 
grees of freedom, and presence of unique corresponding 
gauge-invariant variables) continue to be valid even in 
higher-order perturbations, [11]. 

Whereas, in the synchronous gauge, by imposing (3 = 
in all coordinates (i.e., $ = = (3), from Eq. (48) we have 

C = 0. (53) 

Thus, even after imposing the gauge condition we have 

C Q = £ Q (x), (54) 

which is the remaining gauge mode. Thus, under the 
synchronous gauge, from Eqs. (47), (48) we still have 

a = a — a Q £ a , 5 = S — <5, Q £ Q , k = K — K a £ a , (55) 

where the transformation of 7 is valid to the linear order. 
In this sense variables in the synchronous gauge have re- 
maining gauge modes even after imposing the gauge con- 
dition. 7 has the remaining spatial gauge mode even in 
the linear order, and the other variables have remaining 
gauge modes to the second order. 

C. Transformation between the two gauges 

Using the gauge transformation properties of the vari- 
ables we can translate the equations and solutions in one 
gauge into the ones in another gauge condition. We indi- 
cate the comoving gauge and the synchronous gauge by 
subindices C and S, respectively. To the linear order we 
have 

Pc = is = \x- (56) 
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We present three different ways to reach the transforma- 
tion properties. 

First, we consider a transformation from the syn- 
chronous gauge (unhat) to the comoving gauge (hat). 
From Eq. (47) we have 

Sc = 5s-Ss,aC K C = K S - K S , a ^ a ■ (57) 

We need to determine the gauge transformation function 
£ Q , apparently, only to the linear order. From Eq. (48) 
we have 

£a = IS.a- (58) 

Thus, Eq. (57) becomes 

5c = $S - 6s,aj's, K C = KS - KS,a7s": (59) 

where we used Eq. (56). 

Second, we consider a transformation from the comov- 
ing gauge (unhat) to the synchronous gauge (hat). From 
Eq. (47) we have 

Ss = Sc-S c ,aC k s = k c - K C , a ^ a ■ (60) 
From Eq. (48) we have 

£,a = -7S,a- (61) 

Thus, Eq. (60) leads to the same results in Eq. (59). 

Finally, the gauge invariant combination in Eq. (49) 
provides a simpler derivation. From the gauge invariance 
of combinations in Eq. (49) we directly have Eq. (59). 

Using these gauge transformation properties in Eq. 
(59) we can derive Eqs. (38), (39) from Eqs. (29), (30), 
and vice versa. 



V. DISCUSSION 

In this work we have compared the general relativis- 
ts weakly nonlinear cosmological perturbation equations 
in two different gauge conditions. In our previous works 
we have successfully shown that, except for the coupling 
with gravitational waves, the relativistic perturbation 
equations of a zero-pressure irrotational fluid coincide 
exactly with the Newtonian ones to the second order 
in perturbations. Such a relativistic-Newtonian corre- 
spondence was available in our special comoving gauge 
condition in which all the variables can be equivalently 
regarded as gauge invariant ones. In this work we have 
compared these results with the ones in the synchronous 
gauge. The case in the synchronous gauge was previously 



studied without noticing the similarity or difference of 
the equations with the Newtonian ones to the nonlinear 
orders. 

In this work we compared equations in the synchronous 
gauge with the ones in the comoving gauge and in the 
Newtonian case. Although the variables in this gauge 
have remnant spatial gauge modes due to the incomplete 
gauge fixing of the spatial gauge condition the equations 
are gauge invariant. Ignoring the gravitational waves, 
the equations in the synchronous gauge can be identi- 
fied with the Newtonian hydrodynamic equations in the 
Lagrangian frame to the second order, whereas the equa- 
tions in the comoving gauge can be identified as the New- 
tonian ones in the Eulcrian frame. These Eulerian and 
Lagrangian correspondences can be understood because 
the fluid four-vector in our comoving gauge is in fact 
normal as in Eq. (28) whereas the four-vector in the syn- 
chronous gauge is both normal and comoving (thus La- 
grangian) as in Eq. (37). In our way to clarify the case in 
the synchronous gauge we have addressed and resolved 
several issues related to the two gauge conditions often 
to fully nonlinear orders in perturbations. 
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APPENDICES 

A. Invariance of fluid quantities 

Here, we show that the fluid quantities based on the 
fluid four-vector in Eq. (8) do not depend on the choice 
of g 0a (the spatial gauge condition) to all orders in per- 
turbations. Using a fluid four-vector u a the energy- 
momentum tensor is decomposed into fluid quantities as 
[14,11] 

fab = flUaUb + p (gab + U a U b ) + q a U b + q b U a + TT ab \ (62) 
fl = f ab U a U b , p=^f ab h ab , q a = -f cd U c h d a , 

n ab = f cd h c a h d b -ph abl (63) 

where h ab = g ab + u a u b ; we have u a q a = = u a n ab , 
TT ab = TT bai and 7r° = 0. The variables jl, p, q a and jr a b 
are the energy density, the isotropic pressure (including 
the entropic one), the energy flux and the anisotropic 
pressure (stress) based on the fluid four-vector, respec- 
tively. Let us introduce another four-vector U a with 

U = -a, U a = Q; U° = i U a = -ag v a . (64) 
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Thus, U a is subject to the same conditions as u a in Eq. 
(8), but with possibly different spatial gauge condition 
which could lead to cjif ^ g 0a . The fluid quantities based 
on U a are similarly defined as in Eqs. (62), (63) with U a 
replacing u a : for example, we have jl u = T a bU a U b , etc. 
We can easily show that if p = q a = n a b = we have 

fj u = f ab U a U b = jlu a u b U a U b = fj,u u U U = /i, (65) 

and p u = (}a = — 0) an d vice versa. This result is 
also valid to fully nonlinear order. 

B. Justification of Eq. (36) 

Here, we show that in a zero-pressure irrotational 
medium we can take the original synchronous gauge 
(a = = (3) together with the temporal comoving gauge 
(v = 0) simultaneously to all orders in perturbations. 
This was known in [2], see Sec. 97 in [2]. Here, it is im- 
portant to take (3 = as the spatial synchronous gauge 
although we prefer to take 7 = as the spatial gauge con- 
dition because of the remnant gauge mode in the (3 = 
case. We provide two different proofs based on the ADM 
and the covariant formulations. 

We begin by taking v = and (3 = as the temporal 
and spatial gauge conditions, respectively. In Eq. (7) we 
showed that g 00 = -l/N 2 = -I /a 2 . The spatial gauge 
condition (3 = together with the irrotational condition 
implies gou = N a = 0. Thus, from Eq. (2) of [11] we have 
.9oo = -a 2 (l + 2a) = -N 2 = -a 2 . This implies that we 
have a = 0. 

Now, in the covariant approach, v = and irrotational 
conditions imply ii a — 0. Since u a is the fluid four- 
vector we take the energy frame, q a = 0. The momentum 
conservation equation in Eq. (27) of [11] gives a a = 0. 
The spatial gauge condition (3 = together with the 
irrotational condition implies goa = 0, thus g 0a = as 
well. Since a a = u a - b u b = fg Q = ^g 00 g o, a , d a = 
implies that g o is a function of time only. Thus, we have 
a = 0. 

If we impose a = and (3 = as the gauge condition, 
instead, we have non- vanishing J a or u a , thus nonvan- 
ishing f. In a zero-pressure medium this nonvanishing v 
can be identified as the remnant temporal gauge mode, 
which can be set equal to zero without losing physical 
degree of freedom. 
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